o: 2 



(N 



00 



O 



O 



Compact stars in the QCD phase diagram II (CSQCD II) 

May 20-24, 2009, KIAA at Peking University, Beijing - P. R. China 

http : //vega . bac . pku . edu . cn/rxxu/csqcd . htm 



Vortex Dynamics in Color-Superconducting Quark 
stars: The Re-heating of Magnetars 



Brian Niebergal^''^ , Rachid Ouyed^ , Rodrigo Negreiros^ , Fridolin Weber^ 
^Department of Physics and Astronomy, University of Calgary, 2500 University 
Drive NW, Calgary, Alberta, T2N 1N4, Canada 



Department of Physics, San Diego State University, 5500 Campanile Drive, San 



Diego, California 92182, USA 
Email: bnieber@phas.ucalgary.ca 



rrri. Abstract 

P-1 

hrl \ Compact stars made of quark matter rather than confined hadronic matter, 

_; ' are expected to form a color superconductor. This superconductor ought to be 

Q^i threaded with rotational vortex lines within which the star's interior magnetic 



field is confined. The vortices (and thus magnetic flux) would be expelled from 



H ' the star during stellar spin-down, leading to magnetic reconnection at the sur- 

^ . face of the star and the prolific production of thermal energy. In this Letter, 

we show that this energy release can re-heat quark stars to exceptionally high 
temperatures, such as observed for Soft Gamma Repeaters (SGRs), Anomalous 
^ ■ X-Ray pulsars (AXPs), and X-ray dim isolated neutron stars (XDINs). More- 

^D ! over, our numerical investigations of the temperature evolution, spin-down rate, 

and magnetic field behavior of such superconducting quark stars suggest that 
^^ ■ SGRs, AXPs, and XDINs may be linked ancestrally. Finally, we discuss the 

[^ , possibility of a time delay before the star enters the color superconducting 

^P \ phase, which can be used to estimate the density at which quarks deconfine. 

We find this density to be five times that of nuclear saturation. 



The interiors of compact stars provide a naturally-occuring environment that can 
^ ' be used for studying the properties of ultra-dense baryonic matter. A common means 

j^ ■ of probing this environment is through direct observations of thermal emission from 

the surface of compact stars. By comparing these observations with theoretical models 
one can retrieve key information about the physical processes occuring in matter 
compressed to ultra-high nuclear densities [9l [311 1211 [Ml [321 123] . 

Most attempts to model thermal emission make use of what is called the minimal 
cooling scenario, which involves cooling through the minimum set of particle pro- 
cesses that are necessary to explain the thermal evolution of the majority of compact 
stars. Appropriately, it is used as a benchmark for observations of cooling neutron 
stars. However, there are a number of compact stars (see Table [T|) possessing thermal 
emissions significantly out of agreement with the minimal cooling scenario, indicating 



that other processes may be occuring and need consideration. Some specific classes of 
compact stars that disagree with minimal cooling, and are already distinct for some 
of their other features, are Soft Gamma-ray Repeaters (SGRs), Anomalous X-ray 
Pulsars (AXPs), and X-ray Dim Isolated Neutron stars (XDINs). It is generally ac- 
cepted that SGRs and AXPs are the same type of objects, and it has been speculated 
before that XDINs are also related |29]. Observations indicate that SGRs and AXPs 
are very hot objects. Heating by magnetic field decay in the crust has been suggested 
by Pons et al. [25] as a possible explanation and, at certain times during the star's 
evolution, it has been shown to possibly be the dominant source of heating depending 
on the values assumed for magnetic field decay timescales [Ij. However, even with 
the most liberal decay timescale parameters, crustal magnetic field decay does not 
produce enough heat to explain SGR and AXP observations. In fact, there is no 
physical model that can explain the unusually high temperatures of the objects in 
Table [U Phenomenological studies have been performed [7] , which parametrize the 
magnetic field decay to fit observations, and conclude that some efficient mechanism 
of magnetic fiux expulsion from the star's core is required. Other studies proceed by 
introducing an artifical heat source in an internal layer [12] . 

In this Letter, we expand on the idea that magnetic field decay from compact stars 
causes a re-heating of such objects. As stated above, the net effect of this mechanism 
on the temperature of ordinary neutron stars is much too weak to accommodate 
the temperatures observed for SGRs, AXPs, and XDINs. The situation changes 
dramatically if one assumes that these objects are made of superconducting quark 
matter rather than confined hadronic matter. As demonstrated in our previous work 
[T7] and below, magnetic fiux expulsion from the cores of such stars provides a very 
efficient and robust mechanism that can re-heat compact stars to the temperature 
regime observed for AXPs, SGRs, and XDINs. 

The feature that quark matter forms a color superconductor is critical in our 
study. The condensation pattern that is considered here is the color-fiavor-locked 
(CFL) phase |26l |2l |3l [27] where quarks of all colors and fiavors pair together to 
form Cooper pairs. Henceforth we refer to such stars as color-fiavor-locked quark 
stars (CFLQS) |1|. Because of stellar rotation, CFLQSs develop rotationally-induced 
vortices. The total number of vortices at any given time is A^^, = nQ, with n being the 
vortex circulation and Q the star's spin angular frequency. Hence, as such a star spins- 
down the number of vortices decrease. They do so by being forced radially outward 
[in] and upon reaching the surface are expelled from the star. The magnetic field, 
which exists only inside the vortices, is then also expelled and the subsequent magnetic 
reconnection leads to the production of X-rays [22]. Since the star's spin-down rate 
is proportional to the magnetic field strength, when vortices are expelled, along with 
the magnetic fiux contained within them, the spin-down rate also decreases. Thus, 



^ Other condensation patterns for quark matter fail to produce enough re-heating [SI [10111] while 
pairing and superconductivity in neutron matter is also insufficient [55] , 



Table 1: Compact stars with reported thermal emissions significantly higher than can 
be predicted by the minimal cooling scenario, specifically SGRs, AXPs, and XDINs 
pp. The listed ages are the spin-down ages modified by the vortex expulsion process 
(P/3P), instead of the usually assumed spin-down age (P/2P). 
Name T x 10^ Age 

(K) ( 10=^ years) 

0.15 

2.5 

0.5 

1.3 

6.0 

3.0 

0.73 

153 

47 

4.5 

11.3 

1266 

974 



the spin-down rate and magnetic field strength of a CFLQS become coupled. 
The X-ray luminosity from this vortex expulsion process is given by |16] 
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Lx ~ 2.01 X 10^4 erg s~i r^x.o.i Plu , (1) 

where r]x,o.i is the re-connection efficiency parameter in units of 0.1 and P the spin- 
down rate in units of 10~^^ s s~^. An estimate for the latter as well as the magnetic 
field evolution, both derived from vortex expulsion, is 

where Pq denotes the rotational period of the star at birth, Bq the magnetic field at 
birth, and r is a relaxation time given by 



r = 5xl0Mi^ ^ T^ F^ yrs- (3) 



lO^^Gy /Po\^ / M \ /lOkmV 

One can see from equations (P&E]) that the X-ray emission decreases as the star spins 
down. This X-ray emission is produced on the surface of the CFLQS, which alters 
its thermal evolution. To study this numerically, the general relativistic equations of 



energy balance and thermal energy transport need to be solved. These equations are 
given by (G = c = 1) 



dm 



e.e'' + c.^^\, (4) 

(5) 



IQn'^r^KpJl — 2m /r 

respectively |31j. Here, r is the distance from the center of the star, m,{r) is the mass, 
p{r) is the energy density, T{r,t) is the temperature, l{r,t) is the luminosity, 0(r) is 
the gravitational potential, eu{r,T) is the neutrino emissivity, Cv{r,T) is the specific 
heat, and K{r,T) is the thermal conductivity [31]. The boundary conditions of (jl]) 
and ([5]) are determined by the luminosity at the stellar center and at the surface. 
The luminosity vanishes at the stellar center since there is no heat flux there. At the 
surface, the luminosity is defined by the relationship between the mantle temperature 
and the temperature outside of the star [5] . 

Heating/cooling mechanisms used in our work, other than those included in the 
minimal cooling scenario, are the quark direct Urea process, to emissivity of which is 
on the order of 10^^ erg s~^ cm~^, and the quark modified Urea and Bremsstrahlung, 
which are of order 10^^ erg s~^ cm~^. Due to quark pairing in the CFL state, however, 
the direct Urea process is suppressed by a factor e~^'^ and the modified Urea and 
Bremsstrahlung by a factor e~^^/-^ for T < Tc, where A is the gap parameter for the 
CFL phase and Tc is the critical temperature below which strange matter undergoes 
a phase transition into CFL matter. In this work we have used T^. = 0.4A = 10 MeV. 

The heat produced by vortex expulsion occurs in an emission region just above 
the star's surface. In this region the energy, released by the magnetic field decaying 
after it has been expelled from the star's interior, is deposited. An emissivity for this 
process can be calculated from the energy per unit time (cf. Eq. [1]) per unit volume, 
where the volume of interest is a shell surrounding the star. The width of the shell is 
estimated to be the minimum length of a vortex still inside the star, just before it is 
finally expelled. Such a vortex would be a distance of approximately the inter-vortex 
spacing away from the surface of the star. The shell width is then 

. /P\'/V R \^/V300MeV\'/' 

^ "Jl^i-^ III I I 1 



Afl = 8.42 X lO-km ^-j ^^^j ^— ^J , (6) 

where /i is the average chemical potential throughout the star. 

We have solved equations (jl]) and ([5]) numerically for models whose structure 
(composition and bulk properties) were computed from the Tolman-Oppenheimer- 
Volkoff equations [9l [31]. The results are shown in Fig. [H where the redshifted sur- 
face temperature evolutions for various types of cooling scenarios are plotted. These 



scenarios are CFLQSs (CFL stars with superconductivity and the resulting vortex 
expulsion), and in the shaded region in Fig. [H CFL stars (without vortex expulsion), 
uds stars (strange quark stars without pairing), and neutron stars [31]. The initial, 
or CFLQS birth, spin-periods and magnetic field strengths were chosen such that the 
currently observed SGR/AXP/XDIN values of spin-period, spin-down rate, magnetic 
field strength, and luminosity would be consistent with the values derived from vor- 
tex expulsion (cf. Eq. [2]) . The observed data is taken from Table [H where ages and 
temperatures of SGRs, AXPs, and XDINs are listed. 

Figure [1] shows that CFL quark stars without vortex expulsion cool down too 
rapidly to agree with observed SGR/AXP/XDIN data. One might expect that the 
suppression of emission processes, due to CFL pairing, would keep these stars hotter 
for longer, but pairing also changes the specific heat capacity, resulting in enhanced 
cooling. Standard neutron star cooling processes also lead to stars with temperatures 
much lower (albeit warmer than uds and CFL stars without vortex expulsion) than 
values observed for SGRs/AXPs/XDINs. 

Emission due to vortex expulsion dominates all other processes except during the 
first few minutes. However, we also considered the possibility that a neutron star 
undergoes a phase transition to a CFLQS after a delay. Observational motivations 
for considering a delay are; (i) superluminous supernovae and hypernovae [15]; (ii) the 
large discrepancy between SGR/AXP ages are their progenitor supernova remnants 
[20] . Physically, a delay may be necessary if one considers the following; i) the 
time needed for a newly born compact object to spin-down sufficiently such that the 
center reaches nuclear densities [28], ii) the time for the temperature to reach the 
superconducting critical value (T^) and the resulting vortex lattice to form, iii) the 
nucleation time for strangelets to start fusing together [6j. 

Using the quark-nova model [12] the star is reheated to roughly 10^^ K following 
the transition to CFL matter. This energy is from the release of gravitational energy 
during the collapse as well as the latent heat released when converting from hadronic 
to stable strange quark matter [18] 130] [13] . Within the framework of the quark-nova 
scenario the neutron stars with typical values for spin-period and magnetic field, and 
masses greater than 1.5 solar masses, are most likely to make the transition to quark 
deconfinement and the subsequent color-superconducting phase. Upon this transition 
the magnetic field may be amplified via color ferromagnetism to roughly 10^^ G [TT]. 
We choose initial values for the CFLQS's magnetic field strength accordingly. The 
results are shown in Fig. [2] for CFLQSs born after a delay of tqn = 100 years. 

We have computed the thermal evolution for various types of strange quark stars 
and found that CFL quark stars, possessing a rotational-vortex lattice, are in good 
agreement with observed temperatures of SGRs, AXPs, and XDINs. In our model, 
the CFL star spins-down as a result of magnetic braking and expels vortices, which 
contain magnetic fiux, thus decreasing the magnetic field strength in the magneto- 
sphere, resulting in a lower spin-down rate. Hence, the magnetic field strength and 
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Figure 1: Redshifted surface temperature evolution for CFLQSs (CFL stars with 
superconductivity and the resulting vortex expulsion). The initial spin-periods and 
magnetic field strengths were chosen such that the evolution, due to vortex expul- 
sion, of spin-period, spin-down rate, magnetic field strength, and luminosity would 
be consistent with observations. The shaded region indicates calculations of cooling 
scenarios for various types of neutron stars, uds stars (strange quark stars without 
pairing), and CFL stars (without vortex expulsion). All calculations were done as- 
suming stellar masses of 1.4 Mq and radii of 10.5 km. The observed data is listed in 
Table [D 



spin period become coupled. By including emission from vortex expulsion in our rel- 
ativistic cooling calculations, we have found that the unusually high temperatures of 
SGRs, AXPs, and XDINs can be predicted. 

From previous papers [Ml |T7] we have also shown that by using our model for 
a CFL quark star the evolution of the spin-period, spin-down rate, and magnetic 
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Figure 2: Thermal evolution of a neutron star undergoing a phase transition into a 
CFLQS after a delay. This delay is estimated to be the time needed for the central 
density of a neutron star to reach the critical value at which CFL matter is favoured 
(~ 100 years). Values of the birth temperature for the resultant CFLQS are of 
order 10^^ K, and were estimated using the quark-nova model [T8]. The CFLQS birth 
periods and magnetic fields are as indicated. 



field strength can also be predicted for SGRs, AXPs, and XDINs. The long-term 
evolution of these properties suggests an ancestral linkage between SGRs, AXPs, and 
XDINs. We also note a paper by |il4j, which confirms that the birth statistics of SGRs 
and AXPs are consistent with the number of observed XDINs. Finally, this study 
suggests that a delay time between supernova and quark-nova events (tqn) is possible. 
If the estimated delay of 100 years is correct, then the density at which quark matter 
deconfines can be calculated to be roughly five times the nuclear saturation density 
The quark-nova delay can be measured by considering the discrepancy between 



SGR/AXP ages are their progenitor supernova remnants. Provided one had accurate 
measurements of this age difference, a more precise value for the density at which 
quark matter deconfines could be inferred. 
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